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Abstract 



We deduce the gluino contribution to the three-loop QCD quark mass anoma- 
lous dimension function within the minimal supersymmetric Standard Model 
(MSSM) from its standard QCD expression. This work is a continuation of 
the program of computation of MSSM renormalization group functions. 
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Typeset using REVTe?^ 



The renormalization group method |l|] is a powerful tool for the study of many physi- 
cally interesting quantities within the Standard Model and beyond. Although experimental 
measurements at the highest available energy are consistent with the standard model 0, 
the observed relationship of the strong coupling constant at the Z and the weak angle as 
well as the value of the mb/fnr ratio vis-a-vis the top quark mass remain strong indications 
of a supersymmetric (SUSY) grand unification above lO^^GeV and a SUSY threshold for 
squarks and sleptons in the 0.1 to 1 TeV region. To use the renormalization group method 
to study the above and other quantities one needs to know the renormalization group func- 
tions - (3 - function and quark mass anomalous dimension 7^. In our previous work we 
have calculated the three- loop QCD j3 function with gluino contribution included 0]. In the 
present work we deduce the three-loop quark mass anomalous dimension function including 
the gluino contribution. Here we use the known result for the standard QCD expression 
of the three- loop quark mass anomalous dimension [|]. The two-loop contributions in the 
ratio mfe/m^ turned out to be 20% of the leading contribution For this reason a full 
calculation of the three-loop contribution would be useful. 

The quark mass anomalous dimension is defined as usual 

~ ^ ~di^ ^ 7m(as), (1) 

The renormalization of a quark mass within the MS type framework has the following 
form 



= ZmTn = rn 



(2) 



where B indicates the "bare" mass. The anomalous dimension function is determined by the 
lowest order pole term in the quark mass renormalization constant within the MS framework. 
That is, 

The Oj coefficients are related via a renromalization group equation that serves as a 
powerful check of the calculation. 

/ d \ d 

[imias) - /3ias)-^jai{as) = -a^— 0^+1(0^, (4) 

where the three-loop QCD j3 function including the gluino contribution and ignoring squarks 
is defined as follows 0. 

At^l^ = asP{as), (5) 

where 
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(a) (b) 

FIG. 1. Three-loop graphs giving nontrivial contribution to the quark mass anomalous dimen- 
sion function with the gluino included. Wavy lines denote gluons and the solid loop corresponds 
to quark or gluino. 

with 

A = -yC^A + ^(iV,T+|c^), (7) 
P2 = -^{Cl + f[NfTC^ + +A^NfTC, + (8) 



The MS renormalization of quark mass can be expressed as the following multiplicative 
renormalization 

Zm = X Z^' (10) 

Where Z2 is the quark propagator renormalization constant and Z;^^ renormalizes the 
quark propagator with / %l){x)il){x)dx operator insertion. The gluino contributions to the 
above renormalization constants are in one-to-one correspondence with quark loop graphs 
and differ from them only by color and symmetry factors. Our procedure, as used to de- 
termine the gluino contributions to Z decay at four-loop level and to the /5-function at 
three-loop level, is to decompose the known QCD results into contributions from separate 
graphs. Then one can determine the color factors that relate each graph with an inter- 
nal quark loop to the corresponding graph with a gluino loop. For instance, a subgraph 
consisting a simple quark loop has the color factor 

A^/Tr(T"T^) = NfTd"^ 

while the color factor for a simple gluino loop has the color (and symmetry) factor 



3 



The relative factor 1/2 is due to the Majorana nature of the gluino. Here T"" and F° are 
the gauge group generators in the fundamental and adjoint representations respectively. For 
the gauge group SU(N) they satisfy 



ab 



Tr (T'^T*) = T5' 
Tr(F'^F^) = W'' 

Thus for graphs with a simple fermion loop subgraph one obtains the gluino contribution 
by making the substitution 

NfT ^ NfT + -j-CA 

Here = is the Standard Model limit and Ug = 1 corresponds to the minimal SUSY 
extension with one octet of gluinos. 

The full set of three-loop graphs contributing in the standard QCD Z2 and Z;^^ were 
given in ref. 0. We reanalyse all of the graphs up to and including three-loop level and 
added graphs with gluino propagators. The calculation of SU(N) group factors revealed 
that all but two graphs with the gluino contribution can easily be obtained from the known 
results by simply replacing Nj to Nf + rigCA- However, for the graph of Fig. la we found 
that in order to take the particular topology into account one needs to make the following 
replacement in the standard QCD result for this graph 

27 

The replacement for the graph in Fig. lb looks like 

If the gluino lies above the squark, the current calculation provides the contribution 
from gluinos alone up to three-loop order. This is a gauge invariant subset of the full SUSY 
three-loop graphs and leaves a vastly reduced number of graphs (those with one or more 
internal squark lines) still to be calculated at this order. If the gluino lies lower in mass than 
the squarks, the current calculation produces the full SUSY anomalous dimension of quark 
mass up to and including the three loop order in the region up to the squark mass scale. 

We obtain the following result for three-loop quark mass anomalous dimension function 
with gluino included. 

71 = (13) 

"c^^ + fc7^-y(TiV, + |c^)] (14) 
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73 = 



129 

27 V ^ 2 



108 



27 



(15) 



The eigenvalues of the Casimir operators for the adjoint (A^^ = 8) and the fundamental 
{Np = 3) representations of SUc(3) are 

Ca = 3, Cf = 4/3, and T = 1/2. (16) 

We obtain the following values for the above perturbative coefficients of the 7-function. 

71 = 4 (17) 
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(18) 
(19) 



We see that the gluino gives a substantial contribution to two- and especially three- 
loop levels. Indeed, the gluino contribution reduces the two-loop coefficient by about 10% 
and reduces the three-loop contribution by about 50%. Such a large contribution might 
ultimately be important in phenomenological applications. 

The anomalous dimension of quark mass along with the QCD /3-function determines the 
running of the quark mass. Indeed, at the three-loop level, for the running quark mass one 
has (see, e.g., [|]) 



where, 



0(a,(/i)) 



-2/5: 



mfifJ'i) _ 0(as(/ii)) 
m/(/i2) (f){asifi2))'' 

«s(/i) 

4:71 
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(21) 



The above equation indicates that there will be a substantial shift in the running of a quark 
mass due to the gluino. 

To verify this, we will need the following equations. The running coupling is parametrized 
as follows: 



asjfJ') 1 P2 log L _ 1 

47r ~ (3iL (3fL^ (3fL^ 

,2 \ 



(Pi log' L - Pi log L + PsPi - Pi) + 0(L-^), (22) 



where L = log(/iVA|^^,. 

The general evolution equation for the running coupling to 0{a^) 



has the form 
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FIG. 2. The gluino effect on the running of the b quark mass 
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(23) 



where the superscript n (N) indicates that the corresponding quantity is evaluated for 
n (N) numbers of participating quark flavors. Conventionally, n (N) is specified to be the 
number of quark flavors with mass < /x (< M). However, the eq.(^) is relevant for any 
n < N and arbitrary fi and M, regardless of the conventional specification of the number of 
quark flavors. The log mi/fi terms are due to the "quark threshold" crossing effects and the 
constant coefficients 2/3 = /S^^"^'' — 38/3 = (32^'^^ — (32^'^ represent the contributions of 
the quark loop in the /5-function. The sum runs over N — n quark flavors (e.g., / = 6 if n = 4 
and N = 5). Note that mi is the pole mass of the quark with flavor /. Quark masses can be 
estimated from QCD sum rules. For instance, the b quark pole mass is = 4.72GeV [10 . 

In fig. 2 we show the evolution of mf,(/i) from /i=4.72GeV to Mz- We see that the gluino 
effect is a few percent at Mz which could ultimately be important for grand unification 
studies. If the gluino is not light but is nevertheless below the squark mass, our current 
result will still have a region of relevance. Ultimately, of course, it will be desirable to have 
the full SUSY three-loop effect including squark contributions. 

Analyzing numerical results, we see that the third-loop effect is only a few MeV. Based 
on this, we conclude that the three-loop anomalous dimension with the gluino included is 
a good approximation for phenomenological applications and the error in the perturbative 
evaluation of this quantity can be estimated at few MeV. 

An observable quantity R is invariant under the renormalization group transformations 
and obeys the homogeneous renormalization group equation: 

d ^, , d , ^ d 



/i 



m 



f 



dm 



R{lJ, as, mj) = 



(24) 



The quantity R may denote various cross sections and decay rates that are usually calcu- 
lated using the perturbation theory method. In our previous works we have calculated the 
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gluino contribution to the QCD /^-function at the three- loop level and to the hadronic 
decay rate of the Z boson to the four-loop level Thus, the present work completes the 
evaluation of gluino contributions necessary for 0{ag) renormalization group analysis for the 
above quantity. These results can be used in the renormalization group analysis for other 
quantities, for instance, hadronic decay rates of the r-lepton and various Higgs bosons. 
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